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The properties of the similarity transformation in percolation theory in the complex plane 
Ft | , of the percolation probability are studied. It is shown that the percolation problem on 
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a two-dimensional square lattice reduces to the Mandelbrot transformation, leading to 
^vj a fractal behavior of the percolation probability in the complex plane. The hierarchical 

chains of impedances, reducing to a nonlinear mapping of the impedance space onto 
itself, are studied. An infinite continuation of the procedure leads to a fixed point. It 
is shown that the number of steps required to reach a neighborhood of this point has a 
I , fractal distribution. 
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In percolation theory, an approach based on a hierarchical chain of resistances is used 
to calculate the percolation threshold and the behavior of a system near the percolation 
threshold. The approach is based on a recurrence relation for the percolation probability 
p |]J of a medium or for the resistances of the hierarchically constructed cluster (see, for 
example, Ref. 2). 

A very simple model of a scale transformation for geometric percolation employs a calcu- 
lation of the percolation probability for a d-dimensional cube with edge length 2a and know 
percolation probability a. If one neglects percolation between the diagonally positioned con- 
ducting cubes, the renormalization-group relation for the percolation probability at the nth 

level has the form 
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Q\ . The mapping (1) is a particular case of more general formulations of the renormaliztion- 

group study of percolation transitions. Our objective in the present letter is to investigate 
the mapping (1) in the complex p-plane. This mapping is nonlinear and its general properties 
are described by the theory of fractals [[3] . Another problem is to study a similar mapping 
for hierarchically constructed chains of impedances simulating the high-frequency properties 
q | of a percolation medium. 

Let us consider the two-dimensional (2D) case. The mapping (1) possesses three fixed 
points: 0,1, and p = (a/5 — l)/2. The point p c is unstable: for po > p c , p n asymptotically 
approaches 1, and in the opposite case p n approaches zero. Therefore p c is the percolation 
threshold. The investigation of fractals is based on the study of the very simple nonlinear 
Madelbrot map 

x' = F(x) = x 2 + c. (2) 

It is easy to see that choosing c = — 1 , substituting Eq. (2) into itself, and replacing 
x by — p gives precisely the map (1) with d = 2. Therefore the Sarychev map (1) is a 
particular case of the Mandelbrot map. By studying this map int the complex p-plane we 
can determine the regions from which a conducting (p n — > 1) or nonconducting (p n — > 0) 
medium is reached as n —>■ oo. In terms of the Mandelbrot theory, they are called basins of 
attraction of the corresponding points 
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The attraction basin of the point zero, corresponding to nonpercolation in the 2D case, 
is displayed on Fig. 1. It is distinguished by coloring the plane as a function of the parity of 
number of steps required to reach a circle with radius 310~ 2 around the origin. The basins 
of attraction of other finite fixed points are colored in black. The basin of attraction for zero 
cuts off (0,p c ) of the real axis the part from the segment (0, 1) corresponding to the physical 
region. The percolation system is determined by the segment (p c , 1), which falls in the basin 
of attraction of the point 1 - the black region on the right-hand side. 

The fractal structure of the basin of attraction is formed by regions spanning the points 
that map into zero after many repetitions of the mapping, i.e. , the points are zeros of the 
functions F n (p) = F(F(...F(p))...). This function determines the percolation probability of 
a square with an edge length 2 n . The derivative of F n (p) with respect to p can be regarded 
as the distribution function of the percolation thresholds. F n is a polynomial of degree 2 n 
with even powers. 

It is easy to obtain the asymptotic behavior of F n as Sp = p—p c —* 0. Linearizing the map 
(2) gives Sp' = 2p c Sp i.e. F^{p c ) ~ (2p c ) n . Since 2p c > 1, the derivative grows exponentially 
with n and F n (p) — *■ Q(p — p c ) as n — > oo. 

The general form of this function is, however, quite complicated. It can be shown that 
the general formula for the zeros of this polynomial has the form 



O-l yl + O^V 1 + (T 3 y/.... + (T n , (3) 

where Oi = ±1. Evidently, all zeros before the mappings are also zeros after the mappings. 
The preceding zeros appear with a higher multiplicity. Specifically, if n = 2 k , then the zero 
of the 2 J th mapping (j < k) appears with multiplicity 2 k ~K 

In the limit n — *■ oo the zeros approach the boundary of the basin of attraction, forming 
a so-called Julia set. The fractal nature of the Julia set can be discovered by examining 
the periodic subsequences {o"j}. Each such subsequence approaches its own irrational limit. 
Together the limits all form a fractal boundary. 

It is easy to obtain a physical interpretation of a basin of attraction. Consider the quantity 
Re(ln F n (x)) . This quantity is 2D electrostatic potential of charges which are placed at zeros 
of the function F n (x) and have a magnitude equal to the multiplicity of the zeros. The 
equipotential lines form the basin of attraction. Large, circular regions of influence of the 
charges arise near the largest charges-the zeros of the mappings with small j. As j increases, 
the number of regions increases as the number of zeros and size p n of regions decreases. The 
rate of decrease of the size can be found by equating the potential produced by the charges 
at the points where the boundaries of the corresponding regions touch: 
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whence p n = p$. The choice of the set {u n } equal to 1 gives, in accordance with the 
Lee- Yang theorem, a chain of zeros of a polynomial which approach p c . 

Therefore the percolation probability behaves very simply on the axis in the limit of large 
n, viz. as 0(p — p c ), and in the complex plane it exhibits a complicated, fractal behavior. If 
one intersted only the asymptotic behavior of F n (p) for large n and near a transition, then 
there is no need to know the singularities of this function in the complex plane. However, if 
we wish to describe the function for all n and far from a transition point, then knowing the 
zeros of F n (p) makes it possible to reconstruct the entire function. 

Another example of the fractal behavior is the frequency dependence of the impedance 
of a large, hierarchically constructed circuit in the case of small, real losses. In Ref. 4 
it was found that in a random medium consisting of capacitors, inductors and resistors 
with small resistance, the local electric field is subject to anomalously strong fluctuations. 



For example, in a 2D Dykhne medium consisting of two metals the imaginary pert of the 
effective permittivity at a frequency between their plasma frequencies is nonzero in the 
abscence of local absorption. This causes the local electric field to diverge. We shall show 
that anomalously strong fluctuations in the properties occur even in a regular hierarchical 
chain. Strictly speaking in the abscence of absorption even self-averaging of the system does 
not occur. 

Consider now the electrical circuits in Fig. 2. Let the initial element i be a capacitor 
C with impedance Zi = i/uoC at frequency u and the impedance of the initial element of 
the iteration sequence Zq = —iLu + R, where L is the inductance and R is the resistance. 
Switching to dimensionless variables z n = uCZ n , zq = —iLCu) 2 + RCu, we obtain for the 
circuit in Fig. 2a 

Zn+l = %~ ( — T-. (4) 

{Z„ + I) 

The explicit formula for z n has the form z n = (if n +iZo — f n )/(fnZo + ifn-i), where /„ are 
the Fibonacci numbers 
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The set of both poles and zeros has a limit point — «(\/5 — l)/2. Near this point z n 
converges poorly. 

In the example analyzed above, the final formula reduces to a linear fractional transfor- 
mation. In the general case the map (2) performed n times leads to rational function of 
high order. The result of this mapping in the space Zq can be predicted if the stationary 
points and poles of transformation are known. As a degree of the transformation increase, 
the number of poles increase as 2 n and the convergence of the mapping therefore becomes a 
fractal function of the initial impedance. 

Fig. 2b corresponds to the map 

Zn+l ~ 2z n + i ■ (5) 

Figure 3 displays the impedance after 10 iterations as a function of zq. It is evident from 
the figure that in the abscense of absorption the impedance is a fractal curve. According 
to the definition of zo, the fractal nature is reproduced in the frequency dependence of 
the impedance. The fractal behavior is explained by the fact that random circuits which 
resonate at the frequency uo arise in the system. In the limit of small R, the width of a 
resonance is determined by the quantity R/L. The total number of resonances distributed 
in a characteristic frequency band of the order of y/LC is determined by the number 2™ + 1 of 
the element in the system. As the static conductivity of elements increases, the resonances 
cease to overlap if \ / LC2~ n L/R > 1, and the absorption is smoothed out. 

An important property of a percolation system is the percolation radius r c . The geometric 
percolation radius is determined as the scale on which the probability of two sites belonging 
to the same cluster starts to decrease more rapidly than in power law. The correlation radius 
becomes infinite at percolation threshold. 

In the analysis of the conductivity of a two-phase random system containing phases with 
finite dc conductivities, it assumed that the correlation radius is related to the ratio h = Oxjoi 
of the conductivities: r c ~ h^ 13 . However, the question of to determines the correlation radius 
remains open in the case of an ac field. A necessarily real quantity, which r c is, cannot be an 



analytic function of the impedances. Therefore the function r c {ai,a 2 ) cannot be suggested 
as an alternative. 

In the model corresponding to Fig. 2b, the impedance approaches zero as the number of 
iterations increases. In the limit of a sufficiently large number of steps, the renormalization- 
group equation can be converted into the differential equation dz/dn = —z 2 with the solution 
z = l/(n + A). It is natural to take as the correlation length, is units of n, the value of n 
for which z reaches this asymptotic law. For real z, this occurs for n ~ In z. 

In Fig. 4, the number of iterations at which the impedance drops below 0.01 in modulus is 
displayed as a function of the initial impedance. One can see how this dependence becomes 
fractal as R decreases. 

The correlation radius estimated from the smoothing of the resonances has the form 
n ~ In (l/i?)/m2. The fact that r c becomes infinite can be understood if the system under 
study is interpreted in the limit R — > as a system of coupled, localized states for a photon. 
In such a formulation, the spectrum of the system consists of separate lines, which do not 
vanish in the absence of a relaxational interruption of the phases. At the same time, the 
substitution z — > iz gives the system coupled positive and negative resistances, i.e., the 
system contains amplifying sections, which explains the absence of a thermodynamic limit. 
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Picture signs 

FIG. 1. a-Basin of attraction of the fixed point corresponding to nonpercolation with 
d = 2. The regions where the number of steps required to reach a neighborhood of the point 
is even are colored in black, as are the basins of attraction of other, nonfixed points, except 
infinity. b-Same as Fig. la,d = 3. 

FIG. 2. Circuits of the hierarchical chains a and b corresponding to the expressions (4) 
and (5). 

FIG. 3. Impedance of the chain in Fig . 2b as a function of the imaginary part of the 
initial impedance after 10 iterations. Adjacent curves, except for the bottom curve, are 
raised by 0.2 units relative to one another for better resolution. 

FIG. 4. a-Number of steps required to reach a neighborhood of 310~ 2 as a function of the 
initial impedance for the circuit in Fig. 2b. Adjacent curves, except for the bottom curve, 
are raised by 25 units relative to one another. b-Dual coloring, similar to Fig. 1, of the plane 
of initial impedances. 
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